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Abstract 

Accumulating evidence suggests that cutaneous human papillomavirus (HPV) infection is associated with non-melanoma 
skin cancer (NMSC). Little is known about the natural history of cutaneous HPV. A sub-cohort of 209 men with no NMSC 
history, initially enrolled in the HPV infection in men (HIIVl) study, were followed for a median of 12.6 months. 
Epidemiological data were collected through self-administered questionnaires. Cutaneous HPV DNA was measured in 
normal skin swabs (SS) and eyebrow hairs (EB) for 25 and 16 HPV types in genera p and y, respectively. Any p HPV infection 
was more prevalent in SS (67.3%) compared to EB (56.5%, p = 0.04). Incidence in SS was higher than 20 per 1,000 person- 
months for HPV types 4, 5, 23, 38 and 76. Median duration of persistence of p and y HPV infection was 8.6 and 6.1 months in 
EB, respectively, and 1 1 .3 months and 6.3 months, in SS, respectively. Older age (>44 years vs. 18-30 years) was significantly 
associated with prevalent (SS OR = 3.0, 95% CI = 1.2-7.0) and persistent p HPV infection (EB 0R = 6.1, 95% CI = 2.6-14.1). 
History of blistering sunburn was associated with prevalent (OR = 2.8, 95% CI = 1.3-5.8) and persistent (OR = 2.3, 95% 
CI = 1.2-4.6) p HPV infection in SS. Cutaneous HPV is highly prevalent in men, with age and blistering sunburn being 
significant risk factors for cutaneous p HPV infection. 
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Introduction 

Human papillomavirus (HPV) is a non-enveloped DNA virus, 
primarily infecting stratified epithelium [1,2]. HPV infections are 
broadly classified as mucosal (a genus) or cutaneous (a, P, y, |i, Nv 
genera) [3] . An etiologic role for mucosal HPV infection has been 
established for several cancers, including penile [4] and cervical 
[5,6] as well as a subset of head and neck cancers [7]. In the 
majority of cervical cancers, viral DNA integrates into the host 
genome, further facilitating the transformation of infected cells 
[8,9]. Accumulating evidence suggests that, infection with 
cutaneous HPV types is associated with increased risk of non- 
melanoma skin cancer (NMSC), in both immunocompetent and 
immunocompromised individuals [10-12], perhaps through an 
indirect mechanism. For example, inhibition of ultraviolet 
radiation induced apoptosis has been suggested as a potential 
mechanism of carcinogenic activity of cutaneous HPV [13,14]. 

The earhest evidence for oncogenic potential of P HPV types 
came from studies by Jablonska and colleagues, in 1972, which 
demonstrated the presence of P HPV types in cutaneous lesions of 
patients suffering from Epidermodysplasia Verruciformis (EV) 
[15], a rare hereditary skin disease that often progresses to 



squamous cell carcinoma (SCC) of skin in solar exposed areas. 
Cutaneous HPV DNA has been detected in eyebrow hairs, normal 
skin samples as well as tumor tissues from NMSC cases [16]. 
Several studies have reported a positive association between 
cutaneous HPV DNA and/or seropositivity and NMSC [17-20]. 
Despite its potential role in the development of SCC, very little is 
known about the natural history of cutaneous HPV infection. This 
information is essential to further elucidate the role of cutaneous 
HPV in skin cancer and also to guide future preventive measures 
targeted at reducing the burden of cutaneous HPV infection. 

Few studies have evaluated the prevalence and persistence of 
cutaneous HPV among healthy individuals. Cutaneous HPV is 
thought to be ubiquitous in skin, with hair follicles serving as the 
reservoirs of persistent HPV infection [21,22]. High prevalence of 
cutaneous HPV in skin of infants and young children, indicates 
that exposure to HPV occurs very early in life [23], with 
asymptomatic infection persisting for several years [24]. Preva- 
lence rates ranging from 42% in African population to 70% in 
European and Asian populations have been reported among 
healthy adults [25]. Prevalence of over 90% have been reported 
among immunocompromised patients [22]. 
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Table 1. Baseline characteristics of 209 male participants in the humanpapilloma 


virus infection (HIM) Study, Tampa, Florida. 




Variable 


n (%) 


Age 


18-30 


108 (51.7) 


31-44 


44 (21.0) 


45+ 


57 (27.3) 


Self-ldentifled Race 


White 


155 (74.5) 


Other 


53 (25.5) 


Spanish/Hispanic/Latino 


No 


1 77 (84.7) 


Yes 


32 (15.3) 


Marital status 


Single, Never Married or Divorced/Separated 


148 (70.8) 


Married or Cohabiting, Living Together 


61 (29.2) 


Highest level of education 


Higli school or below 


38 (18.2) 


Vocational school/Some college 


110 (52.6) 


Graduated college/Graduate school 


61 (29.2) 


Skin reaction to season's first sun exposure 


No change in skin color 


37 (17.9) 


Tan with no sunburn 


53 (25.6) 


Mild sunburn that becomes a tan 


80 (38.6) 


Sunburn 


37 (17.9) 


Ever had a blistering sunburn 


No 


103 (49.7) 


Yes 


104 (50.2) 


Lifetime number of blistering sunburns 


None 


103 (49.7) 


1 38 (18.3) 


2 


28 (13.5) 


3+ 


38 (18.3) 


Had an alcoholic beverage in the past month 


No 


35 (16.7) 


Yes 


1 74 (83.3) 


Number of days drank in past month 


0 


35 (19.9) 


1-8 


81 (46.0) 


9+ 


60 (34.1) 


Current smoker 


No 


1 76 (84.2) 


Yes 


33 (15.8) 


Ever smoker 


No 


123 (59.4) 


Yes 


84 (40.6) 


Smoking status 


Never 


123 (59.4) 


Former 


51 (24.6) 


Current 


33 (15.9) 


Ever been diagnosed with an STD 


No 


180 (86.1) 


Yes 


29 (13.9) 
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Table 1. Cont. 





Variable 


n (%) 


Lifetime female vaginal sex partners 


0-1 


42 (20.1) 


2-9 


61 (30.6) 


10+ 


96 (48.2) 


Female vaginal sex partners In past 6 months 


None 


53 (25.6) 


1 105 (50.7) 


2+ 


49 (23.7) 


Lifetime male anal-sex partners 


None 


1 58 (90.3) 


1 + 


17 (9.7) 


Male anal-sex partners In the past 3 months 


None 


182 (98.9) 


1 + 


2 (1.09) 



doi:1 0.1 371 /journal.pone.Ol 04843.t001 



In this study, we describe the natural history of cutaneous HPV 
infection in a cohort of 209 healthy men in Tampa, Florida. We 
evaluated a spectrum of epidemiological measures including risk 
factors associated with prevalence, persistence and incidence 
across two genera of cutaneous HPV infection. We observed that 
cutaneous HPV is highly prevalent in men and that, age and 
blistering sunburn are significant risk factors for cutaneous P HPV 
infection in men. 

Material and Methods 

Study population 

The study population for the current analysis comprises a sub- 
cohort of men who participated in the U.S. site of the HPV 
infection in men (HIM) study, a large, multi-national prospective 
cohort study of the natural history of HPV infection in men 
[26,27]. The HIM study methods have been described in detail 
previously [26,27]. Briefly, between July 2005 and September 
2009, study participants were recruited to the HIM study in 
Tampa, Florida, through mass advertisement targeted to univer- 
sity students, faculty, staff and members of the general population. 
Inclusion criteria were: 1) men aged 18-70 years; 2) reside in 
Florida; 3) have had no prior diagnosis of penile or anal cancers; 4) 
have never been diagnosed with genital and/ or anal warts; 5) have 
not participated in an HPV vaccine study; 6) have no prior 
diagnosis of HIV/ AIDS; 7) have no current penile discharge or 
burning during urination; 8) were not currently being treated for 
an STD; 9) have not been imprisoned or homeless during last 6 
months; and 10) had not been in a drug or alcohol treatment 
program over the last 6 months at enrollment. Between November 
2008 and June 2010, 1,082 participants in the HIM study were 
invited to participate in a supplemental study of the natural history 
of cutaneous HPV infection, requiring additional biospecimen 
collection. In the parent HIM study, the participants were 
followed every six months up to four years of foUow-up. Nine 
hundred and sixty seven men enrolled in the parent HIM study 
had at least one sample of eyebrow hair or skin swab at their 
baseline visit. Of these 967 men, 965 men (99.8%) had three 
samples (one sample of eyebrow hair and one skin swab sample 
each, from sun exposed and unexposed skin), at their baseline visit. 



89.4% and 66.7% men had all three samples at two and three 
visits, respectively. In order to maximize our observation time and 
facilitate estimation of incidence and persistence, the sub-study 
was restricted to those men who had all three samples for at least 4 
visits (n = 211). Of these, 209 men had viable samples and were 
included in the final analysis. 

Ethics statement 

Written informed consent was obtained from all participants in 
the parent study, and a separate addendum to the consent was 
completed by aU participants in the sub-study. The parent study 
and sub-study protocols were approved by institutional review 
boards at each recruiting site, including the Human Subjects 
Committees of the University of South Florida, the Ludwig 
Institute for Cancer Research, Sao Paulo, Brazil, the Centro de 
Referenda e Tratamento de Doencas Sexualmente Transmissiveis 
e AIDS, Brazil, and the National Institute of Public Health of 
Mexico, as described previously [27]. 

Data collection 

Questionnaires. At the enrollment visit, HIM study partic- 
ipants completed a comprehensive self-administered questionnaire 
with information on demographics (age, race, education, and 
marital status), socioeconomic status, medical history, smoking 
status, alcohol consumption, and sexual history. Additional 
questions on risk factors for skin cancer (skin's reaction to season's 
first sun exposure, history of blistering sunburn, etc.) were added to 
the cutaneous HPV sub-study. 

Eyebrow hairs and swabs of normal skin. Although, in 
the parent HIM study, swabs from both sun exposed and 
unexposed skin were collected, only swabs from sun-exposed skin 
were used for HPV DNA measurement in the present sub-study 
due to limited funding. A 5x5 cm area of normal, sun exposed 
skin, from the top of forearm, was pre-wetted (0.9% NaCl) and 
swabbed back and forth five times using a cotton-tipped Dacron 
swab (Digene, Gaithersburg, MD, USA). AU swabs were placed in 
a separate vial and preserved in 500 |al Digene Standard 
Transport Medium (STM) for HPV DNA testing. Three to four 
hairs were plucked from each eyebrow (6-8 eyebrow hairs per 
individual) using disposable tweezers. Hairs with attached follicles 
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Figure 1. Time to incidence of p-HPV infection in normal sl<in 
swabs and eyebrow hairs of men. Kaplan Meier estimate for time to 
incidence of any (3 HPV infection in a sample of 209 men. Participants 
who were negative for all HPV types at baseline were included. Time 
was counted until their first visit with a HPV positive sample or until 
censored. The 'last observation carried forward' approach was taken 
when counting time to incidence (e.g. 0 NA 1 pattern of HPV positivity 
at consecutive visits, the NA was treated as '0' and the time was 
counted in the time to incidence.). This only affected skin swab samples 
since most people were p globin or HPV positive for eyebrows and thus 
their samples were not removed. 
doi:1 0.1 371/journal.pone.01 04843.g001 
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Figure 3. Time to incidence of y-HPV infection in normal sl<in 
swabs and eyebrow Kiairs of men. Kaplan Meier estimate for time to 
incidence of any (3 HPV infection in a sample of 209 men. Participants 
who were negative for all HPV types at baseline were included. Time 
was counted until their first visit with a HPV positive sample or until 
censored. The 'last observation carried forward' approach was taken 
when counting time to incidence (e.g. 0 NA 1 pattern of HPV positivity 
at consecutive visits, the NA was treated as '0' and the time was 
counted in the time to incidence.). 
doi:1 0.1 371/journal.pone.0104843.g003 

were snap frozen in liquid nitrogen. All tissues were stored at 
-70°C until testing for cutaneous HPV DNA. 

DNA extraction from skin swabs and eyebrow 
hair. DNA extraction was performed at the International 
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Figure 2. Time to clearance of p-HPV infection in normal sl<in 
swabs and eyebrow hairs of men. Kaplan Meier estimate for time to 
clearance for all type-specific incident infections a study participant 
(n = 209) had before their last visit with a valid sample. The 'last 
observation carried forward' approach was taken when counting time 
to clearance (e.g. for 0 1 NA 0 pattern of HPV positivity, the NA was 
treated as '1' and the time was counted towards the time to clearance). 
This only affected skin swab samples since most people were p globin 
or HPV positive for eyebrows and thus their samples were not removed. 
dol:l 0.1 371/journal.pone.01 04843.g002 
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Figure 4. Time to clearance of y-HPV infection in normal skin 
swabs and eyebrow hairs of men. Kaplan Meier estimate for time to 
clearance for all type-specific incident infections a study participant 
(n = 209) had before their last visit with a valid sample. The 'last 
observation carried forward' approach was taken when counting time 
to clearance (e.g. for 0 1 NA 0 pattern of HPV positivity, the NA was 
treated as '1' and the time was counted towards the time to clearance). 
doi:10.1371/journal.pone.0104843.g004 
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Table 3. Associations between baseline characteristics and genus (3 HPV prevalence in eyebrow hairs and normal skin swabs of 
men residing in Tampa, Florida. 





Eyebrow 


hairs 




Normal skin swabs 




Variable 


Total 
(n = 209) 


Any |i infection 
(n = 118) n (%) 


Age adjusted 
OR (95 %CI) 


Total 
(n = 156) 


Any p 
infection 
(n = 105) n 

(%) 


Age adjusted 
OR (95 %CI) 


Age* 


18-30 


108 


52 (48.1) 


1.0 (reference) 


72 


44 (61.1) 


1.0 (reference) 


31-44 


44 


27 (61.4) 


1.7 (0.84-3.5) 


33 


19 (57.6) 


0.9 (0.37-2) 


>44 


57 


39 (68.4) 


2.3 (1.2-4.6) 


51 


42 (82.4) 


3.0 (1.2-7.0) 


Race 


White 


155 


87 (56.1) 


1.0 (reference) 


115 


78 (67.8) 


1.0 (reference) 


Other 


53 


30 (56.6) 


1.1 (0.6-2.0) 


40 


26 (65.0) 


0.9 (0.4-2.0) 


Spanish/Hispanic/Latino 


No 


177 


99 (55.9) 


1.0 (reference) 


132 


89 (67.4) 


1.0 (reference) 


Yes 


32 


19 (59.4) 


1 .3 (0.6-2.9) 


24 


1 6 (66.7) 


1 .2 (0.4-3.0) 


Marital status 


Single, Never Married or Divorced/Separated 


148 


82 (55.4) 


1.0 (reference) 


106 


66 (62.3) 


1.0 (reference) 


Married or Cohabiting, Living Together 


61 


36 (59) 


0.9 (0.5-1.7) 


50 


39 (78) 


1 .9 (0.9-4.4) 


Higliest level of education 


High school or below 


38 


26 (68.4) 


1.0 (reference) 


28 


17 (60.7) 


1.0 (reference) 


Vocational school/Some college 


110 


59 (53.6) 


0.6 (0.3-1.4) 


81 


57 (70.4) 


1.8 (0.7-4.5) 


Graduated college/Graduate school 


61 


33 (54.1) 


0.4 (0.9-1.0) 


47 


31 (66.0) 


1.0 (0.4-2.8) 


Skin reaction to season's first sun exposure 


No change in skin color 


37 


20 (54.1) 


1.0 (reference) 


30 


24 (80) 


1.0 (reference) 


Tan with no sunburn 


53 


26 (49.1) 


1 .0 (0.4-2.2) 


41 


28 (68.3) 


0.6 (0.2-1.9) 


Mild sunburn that becomes a tan 


80 


46 (57.5) 


1.2 (0.6-2.7) 


62 


36 (58.1) 


0.4 (0.1-1.0) 


Sunburn 


37 


25 (67.6) 


1 .7 (0.7-4.6) 


21 


17 (81) 


0.9 (0.2-3.9) 


Ever had a blistering sunburn 


No 


103 


54 (52.4) 


1.0 (reference) 


78 


45 (57.7) 


1 .0 (reference) 


Yes 


104 


63 (60.6) 


1 .3 (0.8-2.4) 


76 


60 (78.9) 


2.8 (1.3-5.8) 


Lifetime number of blistering sunburns 


None 


103 


54 (52.4) 


1.0 (reference) 


78 


45 (57.7) 


1.0 (reference) 


1 


38 


24 (63.2) 


1.7 (0.8-3.7) 


27 


21 (77.8) 


3.0 (1.0-8.4) 


2 


28 


16 (57.1) 


1.2(0.5-2.9) 


23 


17 (73.9) 


2.3 (0.8-6.8) 


>2 


38 


23 (60.5) 


1.1 (0.5-2.5) 


26 


22 (84.6) 


3.2 (0.9-10.9) 


Had an alcoholic beverage in the past month 


No 


35 


19 (54.3) 


1.0 (reference) 


27 


1 8 (66.7) 


1.0 (reference) 


Yes 


174 


99 (56.9) 


1.1(0.5-2.4) 


129 


87 (67.4) 


1.2 (0.5-2.9) 


Number of days drank alcohol in past month 


0 


35 


19 (54.3) 


1.0 (reference) 


27 


1 8 (66.7) 


1.0 (reference) 


1-8 


81 


48 (59.3) 


1.3 (0.6-3.0) 


57 


39 (68.4) 


1 .3 (0.4-3.6) 


9+ 


60 


31 (51.7) 


0.9 (0.4-2.1) 


43 


30 (69.8) 


1.4 (0.5-4.2) 


Current smoker 


No 


176 


96 (54.5) 


1.0 (reference) 


131 


89 (67.9) 


1.0 (reference) 


Yes 


33 


22 (66.7) 


1.4 (0.6-3.1) 


25 


1 6 (64) 


0.7(0.3-1.9) 


Ever smoker 


No 


123 


64 (52) 


1.0 (reference) 


88 


59 (67) 


1.0 (reference) 


Yes 


84 


53 (63.1) 


1.2 (0.7-2.3) 


67 


45 (67.2) 


0.7 (0.3-1.6) 


Smoking status 


Never 


123 


64 (52) 


1.0 (reference) 


88 


59 (67) 


1.0 (reference) 



PLOS ONE I www.plosone.org 7 September 2014 | Volume 9 | Issue 9 | e104843 



Natural History of Cutaneous Human Papillomavirus (HPV) Infection 



Table 3. Cont. 







Eyebrow 


hairs 




Normal skin swabs 




Variable 


Total 
(n = 209) 


Any p infection 
(n = 118) n (%) 


Age adjusted 
OR (95 %CI) 


Total 
(n = 156) 


Any p 
infection 
(n = 105) n 

(%) 


Age adjusted 
OR (95 %CI) 


Former 


51 


31 (60.8) 


1.1 (0.6-2.3) 


42 


29 (69) 


0.8 (0.3-1.9) 


Current 


33 


22 (66.7) 


1.5 (0.6-3.4) 


25 


1 6 (64) 


0.7 (0.2-1.9) 


Ever been diagnosed with an STD 


No 


180 


98 (54.4) 


1.0 (reference) 


131 


87 (66.4) 


1.0 (reference) 


Yes 


29 


20 (69) 


1.4 (0.6-3.4) 


25 


1 8 (72) 


0.9 (0.3-2.6) 


Lifetime female vaginal sex partners 


0-1 


42 


23 (54.8) 


1.0 (reference) 


27 


16 (59.3) 


1.0 (reference) 


2-9 


61 


37 (60.7) 


1.1 (0.5-2.5) 


40 


29 (72.5) 


1 .6 (0.5-4.5) 


10+ 


96 


56 (58.3) 


0.6 (0.2-1.4) 


82 


57 (69.5) 


0.9 (0.3-2.7) 


Female vaginal sex partners In past 6 months 


None 


53 


28 (52.8) 


1.0 (reference) 


38 


27 (71.1) 


1.0 (reference) 


1 


105 


65 (61.9) 


1.6 (0.8-3.2) 


81 


63 (77.8) 


1.7 (0.7-4.1) 


2+ 


49 


25 (51) 


1.2 (0.5-2.8) 


35 


1 5 (42.9) 


0.4 (0.1-1.2) 


Lifetime male anal-sex partners 


None 


158 


88 (55.7) 


1.0 (reference) 


118 


81 (68.6) 


1.0 (reference) 


1 + 


17 


12 (70.6) 


1.7 (0.6-5.2) 


14 


10 (71.4) 


1.1 (0.3-4.0) 



OR = odds ratio, CI = confidence interval. *Unadjusted I 
doi:l 0.1 371 /journal.pone.Ol 04843.t003 



gistic regression. All analyses, except for age, were adjusted for age using logistic regression. 



Agency for Research on Cancer in Lyon, France, using the 
Qiagen BioRobot EZ 1 with the EZ I DNA tissue kit according to 
the manufacturer's instructions (Qiagen, Hilden, Germany). 
Briefly, the Dacron swabs were carefully cut into an Eppendorf 
tube using scissors and incubated overnight in proteinase K and 
buffer G2 (Qiagen, Hilden, Germany) at 56°C. An EZl DNA 
Forensic protocol was used to extract the DNA from eyebrow hairs 
according to the manufacturer's instructions. To monitor the 
possible occurrence of cross-contamination between the different 
specimens during DNA extraction, tubes containing buffer only 
were also included (one tube with buffer every ten specimens). 

Viral DNA detection using Multiplex PCR/Luminex 
assay. The detection of viral DNA was performed by a 
multiplex PCR/Luminex assay [28], followed by the HPV typing 
using Luminex beads coupled with (3 or y HPV type-specific 
probes. The assay was performed using approximately 100 ng of 
total DNA and specific primers amplifying a part of the E7 gene 
for 25 HPV types in genus P (5, 8, 9, 12, 14, 15, 17, 19, 20,21,22, 
23, 24, 25, 36, 37, 38, 47, 49, 75, 76, 80, 92, 93, 96) and 16 HPV 
types in genus y (4,48,50, 60, 65, 88, 95, 101, 103,108,109, 
1 1 2, 1 1 6, 1 1 9, 1 2 1 , 1 23) [28,29] . Although, a large number of HPV 
types are being continually discovered, a limited coverage of HPV 
types was offered in the multiplex assay. Two primers for the 
amplification of P-globin were included in the assay to provide a 
positive control for the quality of the template DNA. The positivity 
of the assay was given by the intensity of the fluorescent signal 
detected by the Luminex apparatus and was expressed as the 
median fluorescence intensity (MFI) of at least 1 00 beads per bead 
set. The cutoff was calculated for each HPV-specific probe by 
considering the MFI values obtained with no respective PCR 
product. The cutoff was computed by adding 5 MFI to 1 . 1 X the 
median background value, as described by Schmitt et al. [30] . AU 
MFI values above the cut-off have been considered positive. PCR 



was performed with the QIAGEN multiplex PCR kit according to 
the instructions of the manufacturer. This multiplex PCR protocol 
is highly sensitive, being able to detect only 10 copies of the viral 
genome [28,29]. HPV genotyping was successfully repeated in a 
blind manner, three times in 10 individual subjects, demonstrating 
reliability of multiplex PCR for detection of specific HPV types 
[28]. Similar strategy was used for determining the sensitivity and 
specificity of y HPV assay, showing features similar to those of the 
|3 HPV assay. 

Following PCR amplification, 10 |iL of each reaction was 
analyzed by multiplex genotyping using a Luminex based assay as 
described [3 1] . Results were expressed as the median fluorescence 
intensity (MFI) of at least 100 beads per bead set. Of the 1,385 
samples, 1,234 (89%) were P-globin positive. Thirty (2%) (5-globin 
negative and HPV positive samples were included, while 120 (9%) 
fi-globin negative and HPV negative samples were excluded from 
the analyses. AU the tubes containing buffer only were tested 
negative for HPV DNA and (5 globin. 

Statistical analysis 

Baseline characteristics, including demographics, skin cancer 
risk factors, lifestyle factors and history of sexually transmitted 
diseases, were summarized for the cohort using descriptive 
statistics. HPV prevalence in eyebrow hairs and normal skin 
swabs was estimated by phylogenetic species along with concur- 
rent infections with other HPV types. Type-specific prevalence of 
HPV was defined as the proportion of men who were positive for a 
given HPV out of the baseline cohort with valid samples. Genus- 
and species-specific prevalence was defined as the proportion of 
men who were positive for at least one HPV type in the given 
genus or species. The statistical significance of differences in the 
prevalence of HPV infection between eyebrow hairs and normal 
skin swabs was tested using Fisher's exact test. Incidence (in 



PLOS ONE I www.plosone.org 



8 



September 2014 | Volume 9 | Issue 9 | e104843 



Natural History of Cutaneous Human Papillomavirus (HPV) Infection 



Table 4. Association between baseline characteristics and prevalence of species-specific HPV infection in men residing in Tampa, 
Florida. 







Eyebrow hairs 




Sl<in swabs 






P1 positive n = 63 
Age Adjusted OR 

(95% CD* 


P2 positive n = 83 
Age Adjusted OR 

(95% CD* 


pi positive n = 72 
Age Adjusted OR 
(95% CI) 


P2 positive n = 69 
Age Adjusted OR 
(95% CI) 


Age 


18-30 


1 .0(reference) 


1 .0(reference) 


1 .0(reference) 


1 .0(reference) 


31-44 


1.8 (0.8-4.1) 


1.7 (0.8-3.8) 


1 .0 (0.4-2.4) 


1 (0.4-2.5) 


>44 


2.7 (1.3-6.0) 


2.4 (1.2-5.0) 


3.1 (1.2-7.7) 


3.6 (1.4-9.0) 


Skin reaction to season's first sun exposure 


No change in skin color 






1 .0(reference) 


1 .0(reference) 


Tan with no sunburn 






0.7 (0.2-2.4) 


0.4 (0.1-1.4) 


Mild sunburn that becomes a tan 






0.4(0.1-1.1) 


0.4 (0.1-1.3) 


Sunburn 






1.0 (0.2-4.5) 


0.8 (0.2-3.9) 


Ever liad a blistering sunburn 


No 


1 .0{reference) 


I.O(reference) 


1 .0(reference) 


l.O(reference) 


Yes 


1.3 (0.69-2.6) 


1 .3 (0.7-2.4) 


2.6 (1.2-5.7) 


3.1 (1.4-7.0) 


Lifetime number of blistering sunburns 


None 


1 .0(reference) 


I.O(reference) 


1 .0(reference) 


1 .0(reference) 


1 


1.6 (0.6-4.1) 


1.4 (0.6-3.4) 


2.5 (0.8-7.7) 


2.7 (0.8-8.5) 


2 


1.9 (0.8-5.0) 


1 .0 (0.4-2.6) 


2.6 (0.9-8.0) 


2.3 (0.7-7.4) 


>2 


0.7 (0.2-2.0) 


1.4 (0.6-3.1) 


2.7 (0.7-10.0) 


5.0 (1.4-17.7) 


Female vaginal sex partners in past 6 months 


None 






1 .0(reference) 


1 .0(reference) 


1 






1.9 (0.7-5.1) 


1.6(0.6-4.3) 


2+ 






0.5 (0.1-1.4) 


0.4 (0.1-1.2) 



OR = odds ratio, CI = confidence interval. These factors were selected based on significance in the multinomial logistic regression models. 
doi:1 0.1 371 /journal.pone.Ol 04843.t004 



person-months) was determined for any HPV infection, as well as 
type-specific and species-specific HPV infections in normal skin 
swabs or eyebrow hairs. For analyses of incidence of any HPV 
infection, the participant had to be negative for all HPV types 
(genus P or y) at baseline, and only the first incidence of any HPV 
infection was taken into account. Similarly, for species-specific and 
type-specific HPV infections, the participant had to be negative at 
baseline for that particular HPV species or type, respectively. 
Time to incidence of HPV infection was defined as interval, in 
months, between baseline visit and first visit with a positive HPV 
sample or time point being censored. Due to fewer events, median 
time to incidence of y HPV infection in eyebrow hairs could not be 
estimated. Time to clearance of HPV infection was defined as 
interval (in months) between incident HPV infection and clearance 
of all type-specific HPV infections. 

Persistence of any HPV (genus P or y) or type-specific HPV was 
defined as any or type-specific HPV infection at &2 consecutive 
visits. For a given type of HPV persistence, only participants who 
were positive for the HPV genus or type and had at least one 
follow-up visit after that were included. Kappa coefficient (k) was 
used to determine the concordance of viral infections across 
eyebrow hairs and skin swabs for each P-HPV type. The kappa 
coefficients <0 indicates no agreement, 0-0.20 as slight, 0.21-0.40 
as fair, 0.41-0.60 as moderate, 0.61-0.80 as substantial, and 0.81- 
1 as almost perfect agreement [32]. 

Logistic regression models were used to estimate age-adjusted 
odds ratios (OR) and 95% confidence intervals (CI) for the 



associations between risk factors and prevalence of overall and 
species-specific HPV infections (Pl and P2) (due to sample size 
constraints, risk factors associated with type-specific HPV infec- 
tions were not examined.). Variables that were significandy 
associated with P HPV prevalence in an age-adjusted model were 
further evaluated after stratification by P HPV species. Logistic 
regression was conducted to estimate the association between risk 
factors and incidence as weU as persistence of HPV infection, after 
adjusting for age. Logistic regression analyses were conducted 
separately for HPV infection in eyebrow hairs and normal skin 
swabs. Kaplan-Meier curves were derived for incidence and 
clearance of any HPV infection. An alpha level of 0.05 was 
considered statistically significant. Adjustment for multiple com- 
parisons was not conducted. AH analyses were conducted using 
SAS software, version 9.3 (SAS Institute Inc., Gary, NC, USA) and 
R software, version 2.13.1. 

Results 

The median foUow up time was 12.6 months. As seen in 
Table 1, majority (74.5%) of men self-reported their race as White 
and was between 18 and 30 years of age (51.7%). 

HPV DNA in eyebrow hair and normal skin 

As seen in Table 2, cutaneous HPV prevalence was higher 
among normal skin swab samples than eyebrow hairs, particularly 
for any P HPV type (67.3% in normal skin swabs vs. 56.5% in 
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Table 5. Associations of baseline characteristics with incidence and persistence for any type of p HPV infection in men residing in 
Tampa, Florida. 





Incidence 




Persistence 






Eyebrow hairs Age 
Adjusted OR(95% 
CI) n = 40 


Normal skin swabs 
Age Adjusted OR 
(95% CI) n = 50 


Eyebrow hairs 
Age Adjusted OR 
(95% CI) n = 71 


Normal skin 
swabs Age 
Adjusted 
OR(95% CI) n = 95 


Age* 


18-30 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


31-44 


0.9 (0.3-2.8) 


0.5 (0.1-2.04) 


2.0 (0.8-4.6) 


2.4 (0.9-6.7) 


>44 


0.4 (0.1-1.3) 


0.8 (0.1-4.6) 


6.1 (2.6-14.1) 


1 .2 (0.6-2.4) 


Self-Identified Race 


White 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


Other 


1.3 (0.5-3.3) 


0.3 (0.1-1.1) 


0.58 (0.25-1 .34) 


0.44 (0.2-0.99) 


Spanish/Hispanic/Latino 


No 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


Yes 


0.6 (0.1-2.02) 


2.6 (0.3-23.5) 


1 .04 (0.4-2.7) 


0.8 (0.3-1.9) 


Marital status 


Single, Never Married or Divorced/Separated 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


Married or Cohabiting, Living Together 


1 .6(0.6-4.4) 


0.6 (0.1-2.7) 


1.4 (0.6-3.0) 


0.9 (0.4-1.8) 


Highest level of education 


High school or below 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1 .0 (referent) 


Vocational school/Some college 


4.01 (0.8-20.5) 


3.05 (0.6-14.9) 


1.5 (0.6-3.9) 


1.7 (0.7-4.1) 


Graduated college/Graduate school 


5.9 (1 .03-34.2) 


1 .6 (0.3-7.9) 


1.7 (0.6-5.0) 


1 .2 (0.4-3.4) 


Skin reaction to season's first sun exposure 


No change in skin color 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1 .0 (referent) 


Tan with no sunburn 


0.8 (0.2-2.8) 


3.1 (0.5-19.2) 


0.8 (0.2-2.4) 


1.1 (0.4-3.1) 


Mild sunburn that becomes a tan 


0.8 (0.2-2.5) 


3 (0.6-14.9) 


0.9 (0.3-2.6) 


1.1 (0.4-2.9) 


Sunburn 


0.7 (0.2-3.5) 


1.4 (0.2-12.02) 


0.8 (0.2-2.5) 


1.4 (0.4-4.5) 


Ever had a blistering sunburn 


No 


1 .0 (referent) 


1.0 (referent) 


1 .0 (referent) 


1.0 (referent) 


Yes 


1.1 (0.5-2.7) 


1.58 (0.43-5.89) 


1.45 (0.72-2.93) 


2.27 (1.13-4.54) 


Lifetime number of blistering sunburns 


None 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1 .0 (referent) 


1 


0.4 (0.1-1.7) 


2.4 (0.2-23.3) 


0.6 (0.2-1 .7) 


0.9 (0.4-2.2) 


2 


1.2 (0.3-4.3) 


0.5 (0.1-2.5) 


4.01 (1.3-12.2) 


5.04 (1.3-18.9) 


>2 


2.7 (0.8-9.9) 


Could not be estimated 


1.5 (0.6-4.0) 


5.5 (1.7-18.5) 


Had an alcoholic beverage in the past month 


No 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


Yes 


1.0 (0.3-2.9) 


0.4 (0.04-3.1) 


0.7(0.2-1.8) 


0.6 (0.2-1.5) 


Number of days drank in past month 


0 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1 .0 (referent) 


'1-8 


1.1 (0.3-3.6) 


0.4 (0.04-3.5) 


0.5 (0.2-1.5) 


0.4 (0.1-1.1) 


9+ 


0.7 (0.2-2.6) 


0.3 (0.03-3.9) 


0.9 (0.3-2.9) 


1 .4 (0.5-4.2) 


Current smoker 


No 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1 .0 (referent) 


Yes 


0.5 (0.1-2.05) 


0.8 (0.2-3.7) 


0.9 (0.3-2.4) 


0.4 (0.1-1.02) 


Ever smoker 


No 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


Yes 


0.7 (0.3-1.7) 


0.7 (0.2-2.4) 


1.06 (0.5-2.3) 


1.0 (0.5-2.1) 


Smoking status 


Never 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 



PLOS ONE I www.plosone.org 



10 



September 2014 | Volume 9 | Issue 9 | e104843 



Natural History of Cutaneous Human Papillomavirus (HPV) Infection 



Table 5. Cont. 





Incidence 




Persistence 






Eyebrow hairs Age 
Adjusted OR(95% 
CI) n = 40 


Normal skin swabs 
Age Adjusted OR 
(95% CI) n = 50 


Eyebrow hairs 
Age Adjusted OR 
(95% CI) n = 71 


Normal skin 
swabs Age 
Adjusted 
OR(95% CI) n = 95 


Former 


0.8 (0.3-2.2) 


0.7 (0.2-2.9) 


1.1 (0.5-2.8) 


1 .7 (0.7-4.3) 


Current 


0.5 (0.1-2.02) 


0.7 (0.1-3.5) 


1.0 (0.3-2.7) 


0.5 (0.2-1.3) 


Ever been diagnosed with an STD 


No 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


Yes 


2.1 (0.5-9.3) 


1.3 (0.2-7.5) 


0.9 (0.3-2.4) 


0.5 (0.2-1.5) 


Lifetime female vaginal sex partners 


0-1 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1 .0 (referent) 


2-9 


0.7 (0.2-2.5) 


1.9 (0.3-12.7) 


1.4 (0.5-3.7) 


0.8 (0.3-2.3) 


10+ 


0.7 (0.2-2.2) 


1.1 (0.2-5.2) 


0.9 (0.3-2.6) 


0.8 (0.2-2.4) 


Female vaginal sex partners In past 6 months 


None 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1 


2.1 (0.7-6.2) 


0.1 (0.01-1.3) 


1 .0 (0.4-2.4) 


1 .3 (0.6-3.03) 


2+ 


1.5 (0.4-5.0) 


0.3 (0.03-3.5) 


0.9 (0.3-2.7) 


1 .5 (0.5-4.4) 



OR = odds ratio, CI = confidence interval. 

*Unadjusted logistic regression. All analyses, except for 'age' are adjusted for age. 
doi:1 0.1 371/journal.pone.Ol 04843.t005 



eyebrow hairs). HPV types in genus P, species 1 (5, 12, 21 and 24) 
were the most highly prevalent in normal skin swabs and their 
prevalence in skin swabs was higher compared to that in eyebrow 
hairs (McNemar's p = 0.004). No significant difference was noted 
in the prevalence of (32 HPV types between the two sites of 
infection. Significantly higher prevalence of y HPV infection was 
also observed in normal skin swabs (26.8%) compared to eyebrow 
hairs (15.9%, McNemar's p = 0.002). The concordance between 
HPV DNA in eyebrow hairs and normal skin swabs varied by 
HPV type, with kappa values ranging from no agreement (k = — 
0.022) for HPV type 123 to substantial agreement (k = 0.80) for 
types 47 in pi, 92 in P4 and 112 in y HPV (data not shown). 

As seen in Figure 1 , the median time to incidence of cutaneous 
P-HPV infection was 17.3 months in eyebrow hairs with a slightly 
shorter time to incidence in normal skin swabs (11.2 months). 
Thus, at 12 months, 35% of men had newly acquired cutaneous P- 
HPV infection in eyebrow hairs while 57% men acquired new 
cutaneous P-HPV infection in skin swabs. The median time to 
clearance for any P HPV infection was 6.1 months in eyebrow 
hairs and 6.44 months in normal skin swabs (Figure 2). For y HPV 
infection, the median time to incidence in normal skin swabs was 
13.1 months (Figure 3), while the median time to clearance in 
both, eyebrow hairs and normal skin swabs, was 6. 1 months and 
6.4 months, respectively (Figure 4). 

Incidence and persistence of cutaneous HPV are presented in 
Table 2. In general, incidence and persistence rates were higher in 
normal skin swabs compared to those in eyebrow hairs for P and y 
HPV infection. The top five HPV types with highest incidence of 
infection in eyebrow hairs were 38 (P2), 5 (Pl), 4 (y), 80 (P2) and 23 
(P2), ranging from 5.7 to 7.6 cases per 1000 person months; while 
in normal skin swabs these included 4 (y), 23 (P2), 38 (P2), 5 (Pl) 
and 76(P3), ranging from 21.7 to 28.5 cases per 1000 person 
months. More than 50% of men had persistence of HPV types 25, 
47, 12, 15 and 17 in eyebrow hairs and persistence of HPV types 
20, 21, 47, 23, 50, 108, 25, 37, 76 and 8 in skin swabs. Median 



duration of persistence of P and y HPV infection was 8.6 and 6. 1 
months in eyebrow hairs, respectively and 11.3 months and 6.3 
months, in skin swabs, respectively. 

Risk factors associated with prevalence, incidence and 
persistence of HPV infection 

Men aged >44 years were more than twice and thrice as likely 
to have prevalent baseline P HPV infection in eyebrow hairs and 
skin swabs, respectively, compared to men aged 18-30 years 
(Table 3). A positive history of blistering sunburn was associated 
with a more than two-fold increased prevalence of P HPV 
infection in normal skin swabs, although the association between 
frequency of sunburns and HPV prevalence did not reach 
statistical significance (Table 3). No other demographic, lifestyle 
or sexual history related factors were associated with prevalence of 
HPV. 

Evaluation of species-specific P HPV prevalence indicated that 
men aged >44 years were twice more likely to harbor pi and P2 
types of HPV in eyebrow hairs and more than threefold likely to 
have pi and P2 types of HPV infection in normal skin swabs 
(Table 4), compared to men ages 18-30 years. Similarly, a positive 
history of blistering skin burns was associated with more than two 
fold (OR=2.6, 95% CI =1.2-5.7) and more than threefold 
(OR= 3.1, 95% CI= 1.4-7.0) increased prevalence of pi and P2 
HPV in normal skin swabs, respectively. This association was 
stronger for P2 HPV infection in normal skin swabs (OR =5.0, 
95% CI = 1.4—17.7) among individuals with a history of more than 
2 lifetime episodes of sunburns compared to those who reported 
none. Overall, while the association of these factors with P HPV 
infection did not vary significantly by HPV species, age and history 
of blistering sunburns appeared to have strong positive associations 
with P HPV infection, particularly in normal skin swabs. 

No significant associations were seen between baseline factors 
and incident P HPV infection (Table 5) and prevalence or 
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Table 6. Association of baseline characteristics with prevalence, incidence and persistence of y HPV infection in men residing in 
Tampa, Florida. 







Prevalence 




Incidence 




Persistence 






Eyebrow hairs n = 33 


Normal skin 
swabs n = 56 


Eyebrow hairs 
n = 30 


Normal skin swabs 
n = 90 


Eyebrow 
hairs n = 16 


Normal skin 
swabs n = 49 




Age Adjusted OR 
(95% CI) 


Age Adjusted 
OR(95% CI) 


Unadjusted 
OR(95%CI) 


Unadjusted 
OR(95%CI) 


Unadjusted 
OR(95%CI) 


Unadjusted 
OR(95%CI) 


Age 


18-30 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


1.0 (referent) 


31-44 


0.6 (0.2-1.6) 


0.7 (0.3-1.6) 


0.3 (0.1-1.0) 


0.3 (0.1-0.6) 


2.1 (0.4-11.3) 


0.5 (0.1-1.5) 


>44 


0.7 (0.29-1.7) 


1.1 (0.5-2.2) 


0.6 (0.2-1 .4) 


0.4 (0.2-0.9) 


2.7 (0.7-10.3) 


0.7 (0.3-1.8) 


Self-Identified Race 


White 


1.0 (referent) 


1 .0 (referent) 


1.0 (referent) 


1.0 (referent) 




1.0 (referent) 


Other 


0.6 (0.2-1.6) 


0.6 (0.3-1.2) 


0.8 (0.3-2.01) 


0.7 (0.4-1.5) 




0.4 (0.2-1.2) 


Spanish/Hispanic/Latino 


No 


1.0 (referent) 


1 .0(referent) 


I.O(referent) 


1 .0(referent) 


1 .0{referent) 


1 .0(referent) 


Yes 


0.7 (0.2-2.1) 


0.9 (0.4-2.1) 


1 .4 (0.5-3.8) 


0.8 (0.3-1.9) 


1.2(0.2-7.4) 


0.9 (0.3-3.2) 


Marital status 


Single, Never Married or Divorced/ 
Separated 


1 .0(referent) 


1 .0(referent) 


I.O(referent) 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


Married or Cohabiting, Living Together 


0.5 (0.2-1 .4) 


0.6 (0.3-1.4) 


0.6 (0.2-1.5) 


1.0 (0.5-1.9) 


1.5 (0.4-6.0) 


0.5 (0.2-1.2) 


Highest level of education 


High school or below 


1 .0(referent) 


1 .0(referent) 


I.O(referent) 


I.O(referent) 


I.O(referent) 


1 .0(referent) 


Vocational school/Some college 


0.7 (0.3-1 .8) 


1.5 (0.633.8) 


1 .8 (0.6-5.9) 


1.5 (0.6-3.5) 


1.4 (0.3-6.1) 


1.4 (0.5-4.1) 


Graduated college/Graduate school 


0.4 (0.1-1.4) 


1.2 (0.4-3.1) 


0.8 (0.2-3.1) 


0.5 (0.2-1.2) 


1.3 (0.2-8.4) 


0.8 (0.2-3.05) 


Skin reaction to season's first sun exposure 


No change in skin color 


1 .0(referent) 


1 .0(referent) 


I.O(referent) 


I.O(referent) 


l.O(referent) 


1 .0(referent) 


Tan with no sunburn 


0.9 (0.3-3.2) 


0.6 (0.2-1.7) 


1.7 (0.5-5.5) 


0.9 (0.3-2.6) 


1.0 (0.1-7.4) 


2.5 (0.7-9.0) 


Mild sunburn that becomes a tan 


1.3 (0.4-3.9) 


1 .45 (0.6-3.4) 


1.3 (0.4-4.1) 


0.8 (0.3-2.12) 


0.8 (0.1-5.2) 


2.2 (0.7-7.05) 


Sunburn 


1.5 (0.4-5.3) 


0.73 (0.2-2.1) 


0.2 (0.02-1.7) 


0.3 (0.1-0.9) 


2.5 (0.3-21.4) 


3.2 (0.7-14.1) 


Ever had a blistering sunburn 


No 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


I.O(referent) 


l.O(referent) 


I.O(referent) 


Yes 


0.9 (0.4-1 .9) 


1.1 (0.6-2.0) 


0.6 (0.3-1.3) 


1.0 (0.5-1.9) 


3.2 (0.9-10.8) 


0.8 (0.4-1.9) 


Lifetime number of blistering sunburns 


None 


I.O(referent) 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


1 


0.7 (0.2-2.1) 


1 .5 (0.7-3.3) 


0.2 (0.05-1.1) 


1.2 (0.5-3.1) 


1.7 (0.3-10.8) 


0.5 (0.2-1.4) 


2 


0.6 (0.2-2.2) 


0.9 (0.4-2.5) 


0.5 (0.1-1.9) 


1 .7 (0.6-4.9) 


6.2 (0.8-46.1) 


1.0 (0.3-3.5) 


>2 


1.5 (0.6-4.05) 


0.8 (0.3-2.02) 


1.1 (0.4-3.1) 


0.6 (0.2-1.3) 


3.5 (0.8-15.3) 


1.5 (0.5-4.7) 


Had an alcoholic beverage in the past month 


No 


1 .0(referent) 


1 .0(referent) 


I.O(referent) 


1 .0(referent) 


l.O(referent) 


I.O(referent) 


Yes 


0.6 (0.2-1 .4) 


1 .3 (0.6-3.2) 


0.9 (0.3-2.6) 


0.5 (0.2-1.3) 


0.7 (0.2-2.7) 


1.4 (0.5-4.2) 


Number of days drank in past month 


0 


1 .0(referent) 


1 .0(referent) 


I.O(referent) 


I.O(referent) 


l.O(referent) 


1 .0(referent) 


1-8 


0.5 (0.2-1 .5) 


1.1 (0.4-2.8) 


0.6 (0.2-2.1) 


0.5 (0.2-1.3) 


0.5 (0.1-2.8) 


1.0 (0.3-3.4) 


9+ 


0.9 (0.3-2.4) 


1 .4 (0.5-3.7) 


1 .5 (0.5-4.7) 


0.8 (0.3-2.2) 


0.8 (0.2-3.7) 


2.4 (0.7-8.5) 


Current smoker 


No 


I.O(referent) 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


Yes 


1.1 (0.4-3.2) 


0.7 (0.3-1.85) 


1.1 (0.4-3.1) 


0.9 (0.4-2.2) 


1.2 (0.3-5.7) 


1.5 (0.5-4.7) 


Ever smoker 


No 


I.O(referent) 


l.O(referent) 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


Yes 


1.4 (0.6-3.2) 


0.8 (0.4-1.6) 


1.2 (0.5-2.6) 


0.8 (0.4-1.6) 


3.1 (0.9-10.4) 


1.1 (0.5-2.5) 



Smoking status 
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Table 6. Cont. 







Prevalence 




Incidence 




Persistence 








Normal skin 


Eyebrow hairs 


Normal skin swabs 


Eyebrow 


Normal skin 




Eyebrow hairs n = 33 


swabs n = 56 


n = 30 


n = 90 


hairs n = 16 


swabs n = 49 




Age Adjusted OR 


Age Adjusted 


Unadjusted 


Unadjusted 


Unadjusted 


Unadjusted 




(95% CI) 


OR(95% CI) 


OR(95%CI) 


OR(95%CI) 


OR(95%CI) 


OR(95%CI) 


Never 


l.O(referent) 


1 .0(referent} 


1 .OO(referent) 


1 .00(referent} 


1 .0(referent) 


1 .0(referent) 


Former 


1.48 (0.6-3.7) 


0.86 (0.4-1.9) 


1.2 (0.5-3.1) 


0.8 (0.4-1.7) 


4 (1.01-15.9) 


0.9 (0.3-2.4) 


Current 


1.31 (0.4-4.1) 


0.7 (0.3-1.9) 


1.1 (0.38-3.37) 


0.9 (0.4-2.1) 


2.0 (0.4-10.4) 


1.5 (0.5-4.8) 


Ever been diagnosed with an STD 


No 


1 .0(referent) 


1 .0(referent) 


l.O(referent) 


l.O(referent) 


** 


l.O(referent) 


Yes 


0.6 (0.2-2.4) 


1.05 (0.41-2.6) 


0.2 (0.02-1.3) 


0.7 (0.3-1.9) 




0.2 (0.04-0.9) 


Lifetime female vaginal sex partners 


0-1 


1 .0(referent) 


1 .0(referent) 


1 .0(referent) 


l.O(referent) 


l.O(referent) 


l.O(referent) 


2-9 


1.4 (0.5-3.9) 


0.4 (0.2-1.1) 


0.3 (0.1-1.1) 


0.7 (0.3-1.9) 


1 .82 (0.4-9.3) 


0.5 (0.1-1.6) 


10+ 


0.6 (0.2-2.1) 


1 .0 (0.4-2.4) 


0.6 (0.2-1 .6) 


0.6 (0.2-1.5) 


1.7 (0.3-8.1) 


0.4 (0.1-1.1) 


Female vaginal sex partners in past 6 months 


None 


1 .0(referent) 


1 .0(referent) 


l.O(referent) 


l.O(referent) 


1 .0(referent) 


1.0 (referent) 


1 


0.4 (0.2-1.0) 


1 .4 (0.6-3.0) 


0.8 (0.3-2.1) 


1.1 (0.5-2.3) 


0.3 (0.1-1.1) 


0.7(0.2-1.9) 


2+ 


0.4 (0.2-1 .3) 


0.9 (0.4-2.4) 


1.5 (0.5-4.3) 


1.3 (0.5-3.1) 


0.3 (0.05-1.4) 


0.9 (0.3-2.8) 



OR = odds ratio, CI = confidence interval. 

**No subjects left in the 'Other' racial group and in the 'diagnosed with STD' group. The sample size was too small to conduct age-adjusted analyses of incidence and 

persistence of y HPV infection. 

doi:l 0.1 371 /journal.pone.Ol 04843.t006 



incidence of y HPV infection, in eyebrow liairs or normal skin 
swabs (Table 6). 

Findings on factors associated with persistence of P HPV 
infection were similar to those of prevalent infections described 
above. Age>44 years was associated with a six fold (OR =6.1, 
95%CI = 2.6-14.2) increased persistence of infection with P HPV 
in eyebrow hairs (Table 5) but not in normal skin swabs. History of 
blistering skin burn was associated with persistence of P HPV in 
normal skin swabs (OR= 2.3, 95% CI= 1.2-4.6), with increasing 
frequency of sunburns showing stronger association (Table 5). 
When restricted to White men, history of blistering sunburns was 
associated with 2.7 times increased persistence of any HPV 
infection in normal skin swabs (OR=2.73, 95% CI = 1.19-6.26), 
after adjusting for age (data not shown). As expected, this 
association was stronger than that observed for the overall 
population (OR = 2.3, 95% 01=1.2-4.6) (Table 5). However, 
when analyses were restricted to White men and stratified by HPV 
species, there were no significant associations between history of 
blistering sunburns and persistence of Pl (OR = 0.87, 95% 
CI = 0.29-2.62) or P2 (OR= 1.33, 95% CI = 0.54-3.27) infection 
in normal skin swabs, after adjusting for age (data not shown). 

No significant associations were observed between baseline 
characteristics and persistence of y HPV except a significantly 
increased persistence (OR = 4.0, 95% CI= 1.0-15.9) in eyebrow 
hairs among former smokers compared to never smokers 
(Table 6). 

Discussion 

In this HIM study sub-cohort of 209 healthy men, we observed 
a high prevalence of cutaneous P HPV and y HPV in normal skin 
swabs, higher than that in eyebrow hair samples. Although 
cutaneous HPV seroprevalence rates have been reported previ- 



ously [33-35], seroreactivity is not a direct estimate of HPV 
infection. Hence, we attempted to directly measure type and 
species-specific HPV DNA in normal skin swabs and eyebrow 
hairs. 

Results presented here contrast to that of other studies which 
reported a higher prevalence of cutaneous HPV in eyebrow hairs 
among individuals without a history of skin cancer [36,37] and 
lower prevalence (13-37%) in skin [37,38]. Tarmorshuizen et al. 
reported a prevalence of 54% in eyebrow hairs among healthy 
controls [39], which was similar to our findings. Differences across 
studies could be a reflection of variation in number of HPV types 
examined. For example, while we evaluated 41 HPV types (25 
HPV types from genus P and 16 HPV types from genus y), 
previous studies have evaluated fewer HPV types within each 
genus (6 to 28 HPV types overall across studies) [36-39]. Apart 
from this, methodological differences or even cUlferences in 
geographical location of study populations may explain the 
inconsistencies in HPV prevalence across studies. Indeed, HPV 
prevalence rate has been shown to vary by country [25,36]. In our 
study, while we observed a strong association between history of 
blistering sunburn and both prevalence and persistence of P HPV, 
no association was seen between tanning ability and P HPV 
infection, indicating that exposure to ultraviolet (UV) radiation 
rather than host genetics may be more important determinants of 
cutaneous HPV infection. Thus, geographical variation in UV 
radiation might have contributed to variable prevalence estimates 
in the literature. While we were unable to measure HPV 
prevalence in sun-unexposed skin, our previous study of NMSC 
cases indicates that the prevalence of P-HPV is higher in sun- 
exposed skin swabs (100% prevalence) compared to sun- 
unexposed swabs (95% prevalence) among NMSC cases [16]. 

Older age and history of blistering sunburns were significandy 
associated with both prevalence and persistence of HPV but not 
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with incidence of cutaneous HPV infection. Our findings are 
consistent with previous reports on the association between age 
and prevalence of HPV [36,39] . A positive history of Mistering 
sunburn was associated with a more than two fold increased 
prevalence of (3 HPV infection in skin, with a stronger association 
for P2 HPV infection. While an inverse association between 
frequency of sunburns and P HPV seropositivity has been reported 
previously [33], we found a positive association between lifetime 
history of one blistering sunburn and prevalent f> HPV infection in 
normal skin swabs. We also observed a strong association between 
frequency of blistering sunburns and persistence of P HPV 
infection, likely explaining the elevated prevalence observed here. 
UV exposure mediated immunosuppression may predispose 
individuals to a higher risk of cutaneous HPV persistence and 
hence prevalence [40]. Indeed, higher prevalence of cutaneous 
HPV has been observed among individuals who reported working 
outdoors for longer duration [41]. 

The findings of our study should be interpreted with caution. 
The external validity of these findings is limited to men. Future 
studies should include women, especially given the previously 
reported differences in seroprevalence of specific types of 
cutaneous HPV by gender [35]. Our findings are generalizable 
to non-Hispanic Whites only. However, our study population of 
Florida residents is at higher risk of NMSC compared to other 
U.S. states, corresponding to Florida's higher UV radiation index 
[42] , and hence may represent a population that could be targeted 
for novel NMSC prevention measures. To maximize the sample 
size, persistence of HPV infection was determined using both 
incident and prevalent HPV infections. Therefore, the true 
duration of HPV infection is unknown, given the prevalent cases 
at baseline. Finally, we examined only 25 HPV types in genus P 
and 16 HPV types in genus y. Given the large number of HPV 
types that are being continually discovered, our results on overall P 
and Y HPV infection may not reflect the true associations that can 
be affected by unexamined and unknown HPV types. 
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